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ABSTRACT 


Extensions of the investigations aimed at preventing aircraft destruc- 
tion through the effects of hydraulic ram were conducted. By analyzing 
spark shadowgraphs the existence of a family of transient shock waves 
was confirmed. These shock waves were created by a surface anomaly in 
a projectile and are not unlike those found in fragmentary warheads. In an 
attempt to measure wall deflections subsequent to projectile impact, a 
system based on optical fiber proximity probes was investigated. Mathe- 
matical models based on grid distortions were used as a means of calculat- 
ing shock wave pressures. The AFIT thesis correlating cavity size, energy 
and radius to underwater explosions was verified from an analysis of 16 mm 


films. 
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I, INTRODUCTION 


The term hydraulic ram has been historically associated with the 
various waves causing, usually, catastrophic failure of an aircraft fuel 
cell when penetrated by a projectile. 

When a projectile or fragment with high velocity impacts the surface 
of a liquid medium, a disturbance is produced which merece throughout 
the medium. The disturbance is characterized by a region of increased 
pressure and density together with a small forward motion of the fluid. 
When the energy of the disturbance is large, the disturbance travels as a 
surface of discontinuity and is called a shock wave. (Ref. 1). Typically 
the shock wave phenomenon has been identified as the first phase of hy- 
draulic ram. 

As projectile penetration continues, a pressure field is created that 
tends to-displace the fluid from the projectile, transforming projectile ~ 
~ energy to kinetic energy of the fuel. This is the beginning of the cavity 
phase. | 

The phenomenon associated with these two phases has been studied 
for as long as three wa i.e. 1947 through present. Recent investigations 
into fragment-induced shock waves and cavities have beeninitiated, the 
results of which are leading to a more clear understanding of hydraulic 


ram. 





II, SHOCK PHASE EXPERIMENTS 


A. WAVE GENERATION IN SHOCK PHASE BY FRAGMENTARY PROJECTILES 

Previous experiments conducted at NPSinvolving the impacting of 
0.22 caliber projectiles into a water-filled tank produced shadowgraphs 
such as those shown as Figures land 2, Figure 1 is the projectile in free 
flight , whereas Figure 2 shows the family of waves created afterimpact. 
The flow features visibié in Figure 2 are considered to be: 

1. projectile wake, 

2. reflected bandwave, 

3. reflected bow wave, 

4. incident wake recompression shock, 

5. incident bandwave, 

6. transient entry shock, 

7. secondary transient shock wave (see discussion below). 
In Figure 2, items 6 to 9 occur in the water internal to the tank and items 
1 to 5 occur in the air external to the tank. The dense, thick region in 
each shadowgraph is the test tank wall. Region 8 is the projectile's 
cavity and item 9 is the grid distortion used to determine the medium's 
density. The spark shadowgraph method was used as it has the advantage 
of showing the projectile-induced waves at one observation. Of particular 
interest was the formation of the wave 7, in Figure 2, which was investi- 


gated as follows. 
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A randomly selected, factory loaded 0.22 cal projectile has the geom- 
etry shown in Figure 3. The volume of the projectile's nose was calculated 
using iccestoienran cated cones while the portion aft of the shoulder was 
calculated assuming a solid cylinder. Assuming a constant velocity rane 
ately after penetration into the tank, the rate of volume penetration was 
plotted against distance along the projectile's centerline, Figure 4. The 
discontinuity in Figure 4 appears at a longitudinal distance of 0.20 inches 
and corresponds to the distance where the projectile itself has a discontin- 
uity, i.e., the shoulder. This shoulder is the origin of bandwave shown 
in Figure 1. This suggests that although the Surface irregularity is small 
(0.020 inches), it is still capable of creating secondary transient shock 
waves. It is concluded that the wave 7 in question on Figure 2 is indeed 
a secondary transient entry shock. Reference 2 describes the deformation 
of several previous test rounds. Considering their fragmentary shape it 


could be expected that fuel cells punctured by fragmentation warheads 


would be subjected to a myriad of such secondary shocks. 


B. PRESSURE BY GRID DISTORTION 

Figure 5 is a shadowgraph of a fine-mesh grid taken while the fluid 
was atrest. No grid distortions are apparent. The projectile's path is 
from the left to right through the center of the shadowgraph. After penetra- 
tion the grid in this area is distorted. This region was used in the subse- 
quent analysis. 

As the projectile penetrates the water it displaces fluid by compressing 


it to very high pressures. This disturbance is the transient shock 6 in 
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Figure 2. As penetration continues, the transient shock continues to receive 
energy until the projectile becomes fully wetted; this occurs at 23 micro- 
seconds in Figure 4. A rarefaction wave then forms that tends to dissipate 
the shock wave. (Ref. 3). The localized high shock pressures are known 
to be of short duration, usually tens of microseconds, and are completely 
dissipated by the time significant projectile penetration occurs. 

The spark shadow-graph could be set off with any time delay up to 
approximately 30 microseconds. ‘Each shock wave in the shadowgraph re-~ 
veals a distinct presence-of-light band and anabsence-~of-light band. The 
concentration of light around the projectile's transient entry shock is best 
explained by considering the pressure and density increase due to the shock 
wave. Similar optical distortion at the nose of the fragment is due to high 
density at the nose. See Figure 2. The OR is known to reach values 
in excess of 4000 kg/cm (56,890 psi), which is sufficient to produce a 
change in the medium's refractive index. The width of the dark band varies 
with the pressure in the wave. The fractional change in refractive index, 
€ , has been studied by Poindexter and Rosen and was determined to vary 
as €&€ =1.06x 107°p where p is pressure in kg/cm (Ref. 1). This conforms 
to an isothermal change, whereas a change produced bythe transient shock 
is more nearly adiabatic (Ref. l). A value of E€= 107°p has been justified by 
McMillen and was used throughout this analysis (Ref. 1). It was also as- 
sumed that in a homogeneous medium such as water, light rays have the 


form of straight lines. 
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Seven shock pressures were extracted from Figure 6. When the pro- 
jectile had penetrated three cm into the tank the pressure was noted as 
2700 ka/em*. Define ©, and @, as the angles which the incident ray and 


refracted ray make with a normal to the shock wave surface (Fig. 7). 





Fig. 7 Light rays in homogeneous regions 
bounded by a shock wave. 


By geometry 8, is given by (1) 
sin 9,= © 


R 
@,= SIT, (| 


Using Snell's law the angle 9, can be introduced 
a 
ny Sin 07, = MSHS, 
sin 87 = (sun 6, 
Ny 


~~ 


From the previous paragraph the change in refractive index can be expressed 
as €= ATL = Uo hy 
TN 

Ny 
on (© pe hy (\+e) 


From equations (1) and (2) 


*nei-(tKE) 


". Q)= sui! (ie | (3) 
and ZAG@= (Or 1) + (8-8, = Be 8, \ (4) 





(2) 
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Note that the changes in angle entering and leaving the shock wave are 
additive. If the density is homogeneous within the region bounded by the 
shock wave, as is assumed here, then the high pressure region behaves 
as a spherical lens. Forge = 0.90 from equation (]l), .-F is found to be 
5.7409. For a pressure of 2700 kg/cm2, ®, is calculated to be 5.588° by 
equation (3); consequently ZAS = 0.303°. Even for high pressures the 
angular deflections are small. Calculated values for the remaining six 
pressures are presented in Table 1. These tabulated results appear reason- 
able. 

The grid distortion analysis assumes rays have spherical symmetry 
i.e. the refractive index depends only on the distance r from a fixed point 


0; see Figure 8. 





Fig. 8 [Illustrating Bouguer's formula nd=constant 


For rays Situated in a plane through the origin, nrsinf=c, where c is 
a constant. With d=rsing, nd=c 
This is known as the formula of Bouguer. From geometry, if (r,@) are 
the polar coordinates of a plane curve, the angle db between the radius 


vector to a point P on the Same curve and a tangent at P is given by 


suc = r(e) 


| r’(6) ie (S% y 
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Index of Refraction 
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for various shock wave pressures 


Exit angle, degrees 


Nondimensional radius 
Entry angle, radians 
Exit angle, radians 
Entry angle, degrees 
2(ANS1-ANS2) 


Table 1 
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Angular deflection of light ray by a spherical shock wave 


for various shock wave pressures 


Nondimensional radius 
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for various shock wave pressures 


Nondimensional radius 
Entry angle, radians 
Exit angle, radians 
Entry angle, degrees 
Exit angle, degrees 
2(ANSI-ANS2) 


Table 1 
RR 
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Thetaz2 
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It is possible to eliminate by introducing the constant c, as can be seen 


by the following manipulations: 
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The assumption of constant pressure behind the shock wave is not 
realistic. At this point, the assumption of uniform conditions behind the 
shock wave can be abandoned by introduction of equation (5). In equation 
(5), c is determined by the point where the ray intersects the shock wave. 
The index of refraction, n, which is related to pressure by E= 107°p, can 
be made a function of the radius. 

To adapt equation (5) to the specific case of the hydraulic ram shock 
phase itis aeen Stee to carefully define the various geometric parameters 
involved. Figure 9 illustrates the angles and radii necessary to apply 
equation (5). Note that the diagram is symmetric about r* at angle a. 

A functional form for n(r) can be assumed. For example, a parabolic 
relation would be 

n(r) = ar¢ + br+c 
In order to evaluate the constants a, b, and c, the integral of equation. (5) 
must be-evaluated and iterated with geometric relations. This is quite, 
complicated . In lieu of a parabolic relation a linear formula is used='-®his '->" 
has the merit of a straightforward solution. One of the two parameters, i.é. 
A, can be determined leaving B to be found as the following manipulations 


show. 


Assume (Fig. 9) 


N= A+BL re yo, Re Ny 
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K= n+ i ae 
t bt y! ') 
From the law of Bouguer 


KV SUA = constant 


Kivsug, = nty* such *= uty* 
eee tm 
nv SUG ,= |e B(x - Nal 


From analytic geometry and Snell's law 


SUA, = Wr 
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From Figure 9 with X14 ae WM 7 known 


Yot 4 ty = eis h, 


Yo FY e"su(6- a) = r/ Su, (Bi- *) 
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$i 7+ b= T+ 40 
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ae Low eps») Fon a 


Uy ltr r' tos (@ X)| rane 
and 


Uy =y! Sm(O-s) ale rv! 05 (B,-2)) hom A (17) 
With the following substitutions, the solution for © is 


b= “| ay 
/ , L 
| (ae ZABr4Br’) rr — | t 


C= ror’ Sug, 


B= B 
" ( 
B= a a 
y=r% C[ BY 2AB y+ Kr eV 
la. { 
wWyK od, =|. nyt O( “ir (18) 
} te 
ao, “\ ee ee ee ee (19) 
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Using a pressure behind the shock wave of 400 kg/cm , Ng was 
found to be 1.335. See Figure 6 for data used in these calculations: 
r' is 7.62 cm. As shown in the Appendix the value for B is -0.04947. 
Using the equation for the linear pressure distribution it is calculated 
that pis 3646 kg/cm4 at ~ 0.1267. The pressure increases toward 
the projectile nose. This result is consistent with Figure 29 of Holm's 
thesis (Ref. 2). 

Results from fluid mechanics indicate that linear relations for p(r) 
and n(r) are not appropriate. Consequently the calculated pressure at 
Ty, vlhe 0.1267 is probably not correct. This analysis gives the correct 
trends. The difference between experimental and calculated values is 


considered to result from the assumed distribution for n. 


C. WALL RESPONSE BY FIBER OPTICS PROXIMITY PROBES 

‘The structural response of the entrance wall during the shock phase 
has been investigated and is presented in Refs. 4,5, and 6. Since the shock 
wave is a pressure disturbance propagating through the fluid medium, it is 
also by definition a stress wave. Types of stress waves include compression, 
tension and shear waves. The compression and tension waves involve posi- 
tive and negative pressures, respectively, where the media particle motion 
is in the direction of wave propagation. The shear wave involves particle 
motion that is normal to the wave direction. Compression and tension waves 
are also known as dilatation waves while shear waves are sometimes called 
distortion waves. Tension waves that reduce compressed areas are also 


known as rarefaction waves (Ref. 3). Loads transmitted through the tank's 





structural members are initially transmitted as stress waves of both dilata- 
tion and shear form. The outward motion ofthe entrance wall somewhat 
relieves the pressure disturbances (stress waves) acting on the wall. To 
measure the outward displacement of the wall a system based on optical 
glass fiber bundles was constructed. The optic fiber bundles were arranged 
coaxially with the transmitting bundle on theinside. This light was then 
directed onto the aluminum entrance wall. The reflected light was collected 
and transmitted by the output fiber bundle to an RCA 931A photomultiplier 
tube and a Pacific Photometric Model 1] laboratory photometer. The fiber 
bundles were three feet in length and were comprised of three-mil flint 
glass fibers with an index of refraction of 1.62 anda cladding glass of in- 
dex 1.52. The corresponding numerical aperture was 0.56 with an accept- 
ance conehalf angle of approximately 34 degrees. The glass bundles were 


secured in a mounting bracket and placed over the entrance wall on a radial 


line emanating from the projectile impact point. The mounting bracket al- 


“lowed-for‘positioning the bundles up-to a maximum of. nine inches-from the.,,,...., u, 


‘ae? ath 


impact point. To increase the intensity of reflected light the entrance wall 
was hand polished with #500 Wetordry Tri-M-ite paper, then with Crystalbay 
Crocus cloth. The signal Output of the photomultiplier was transmitted to a 
dual beam oscilloscope fitted with a Polaroid camera. . 

Initial calibration of the system was accomplished by mounting the 
fiber bundles in a platform micrometer. The light was reflected from a 
circular aluminum disk of the type uSed for the entrance wall. This disk 


was Similarly polished. A calibration curve of volts versus the distance 


separating the disk and glass bundles was obtained, Figure 10. 
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_ Subsequent investigation of this system, modified to include a 
"chopper wheel", revealed an unacceptably low rise time for the photo- 
multiplier tube. A Pacific Photometric Model 62 multiplier tube (RCA model 
6199) was substituted; a Spectrum Physics model 132 HeNe laser replaced 
the light bulb. The inherent noisein the lattersystem obliterated the os- 
cilloscope trace, rendering this optical fiber concept useless without 


further development to improve signal-to-noise ratio. 
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Fig. 10 Calibration Curve for Optical Fiber Proximity Probe 
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III. CORRELATION OF OSCILLATION PERIOD, 


CAVITY SIZE, AND ENERGY DISSIPATION DURING CAVITY PHASE 


A. TECHNICAL APPROACH TO CORRELATION 

In displacing water from its path, the projectile imparts a radial 
velocity to the water. The momentum associated with this velocity field 
is sufficiently high to cause the flow to separate from the projectile and 
to continue to expand causing a cavity. Fluid pressures caused by the 
tank boundaries eventually halt the cavity's expansion. At this point 
most of the kinetic energy initially associated with the velocity field has 
been converted to stored energy, gravitational potential energy and de- 
formation of the tank's wall. These forces cause the cavity to start to 
collapse radially and, as the cavity shrinks, the fluid velocity increases 
resulting in the conversion of stored energy into kinetic energy and work 
in compressing gas trapped in the cavity. Since projectileimpact was near - ---:- 
the center of the tank the cavity would first neck down until anearly spheri- 
cal shape resulted. The sphere would then tend to collapse to a point until 
the pressures in the cavity halted this inward motion. Air which entered 
the cavity through the entrance wall hole is trapped and compressed to a 
high pressure at the instant of minimum cavity volume. A high pressure 
shock wave then propagates radially outward fromthe collapse point. After 
the pressure wave generated by the collapse subsides, the compressed air 
causes the cavity to re-expand. This cyclerepeats several times until all 
the stored energy is dissipated (Refs. 6 and 7). Various authors, Lundstrom, 
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Stepka and Cardea, et al., have likened the cavity phase to underwater 
explosions. Using the method presented by Cardea and Cole, an analysis 
was made of a high speed film showing 0.30 cal and 14.5 mm rounds im- 
pacting a 6061T6 aluminum entrance wall. The film was obtained from 
Naval Weapons Center. These values were then compared with the analyti- 
cal solutions obtained by Holm's research at NPS. Two significant varia- 


tions were noted. Reference 7, equation (7), presents the maximum cavity 
Py Lo 
Po 


accurate. ‘his difference accounts for an 8% variation in tabulated results 


radius as O.< S9T( A constant of -244 is considered to be more 
shown as Table 2. Secondly, this paperconsiders the effects of projectile 
impact through a pre-punched aluminum entrance wall whereas the 0.30 
cal and 14.5 mm rounds impacted a solid tank wall. For the latter the 
shock wave is generated by the impact of the projectile on the wall vice 
the fluid medium. The generation and propagation will be identical ex- 


cept for wall density and sonic velocity effects. This problem is addressed 


». in Reference. 3. The analysis of the 14.5 mm round must be considered as... 


approximate as the maximum cavity radius to be measured extended beyond 


the frame size of 16 mm film. 


B. DATA OBTAINED FROM MOVIE FILM 

For each sequence of the 16 mm film obtained from the Naval Weapons 
Center, China Lake, the maximum cavity radius and time for one complete 
oscillation were measured. Using underwater explosion theory (Ref. 8) 


the energy required to expand the cavity was obtained first as a function 
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of cavity radius and second as a function of the cavity period. A theoret- 
ical value for the cavity radius was also calculated. These results were 


then compared with those obtained by Cardea (Ref. 7). 


C. RESULTS 

Table 2 lists the results obtained from the above measurements. The 
conclusion equating the cavity phase of hydraulicram to underwater ex- 
plosions appears valid. 

For a comparison of the time required for one complete cavity oscilla- 


tion, Reference 2 presents 


2 Poo \+ 
Exal 


where Ory,is the meaSured cavity radius. FromReference 8 this time is 


given as 


\ + 
: Py “e 
where 


Po= 2210 psf 
Ory 1.04 slugs 


Evaluating the constants yields 


4{= 3.L0y (Fe)* 


Po\'1 
£= W84Qm [42 

Po 
respectively. It can be seen that the only difference in characteristic 


times is in the constants. Based on results ofanalysis of the film, 


equation is considered more correct. 


33 





0.30 cal 


Shot number 


nn Pw NHN 


14.5 mm 
Shot number 


Oa fh WD 


e 


millisec's 


24.8 
Dbbag | 
2 oa 
28.8 
Tal ealk 


46.9 
48.3 
47.4 
44,4 
45.9 
43.9 


am 


measured 


Ft 


~485 
2000 
632 
~656 
632 


Zs 
1.094 
PSO Oe 
1.234 
1,266 
1,234 


am 


calc 
Ft 


2455 
404 
.ol6 
2029 
034 


. 861 
. 887 
.870 
yells 
836 
. 806 


Ec=f(Am) 
Ft-lbs 


874.78 
610.68 
IZ7 2.9 
370.0 
1413.3 


hohe lee 
6462.3 
ol07 27 
SUL Tg 
5402.3 
4852.2 


Ec=f(T) 
Ft-lbs 


874.8] 
610.70 
127 2a 
Si Or! 
1413.3 


Sellars i 
6462.5 
6107.7, 
5020.1 
9402.5 
4852.4 


Table 2 Cavity Period, Maximum Say Radius, Cavity Energy 


from Movie Film 


0.50 cal 


Shot Number millisec's 


wn BB W bd 


Table 3 Cavity Period, Maximum Cavity Radius, 


a 


24.8 
26.2 
OS, 
BOOS, 


am - 


observed 


e ® 
wnaoa mo 


from Cardea (Ref. 7) 
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am 


calc 


-494 
aA 
2936 
936 


EC 
Ft-lbs 


S200 
1032.0 
rt 20 
0 


Cavity Energy 





VI. CONCLUSIONS AND PROPOSED RESEARCH 


A. CONCLUSIONS 

Secondary shocks will occur in fuel cells penetrated by fragmenta- 
tion warheads. 

The mathematical model proposed for analyzing the optical grid dis- 
conten yields correct trends. Further development of an assumed diSstri- 
bution for the index of refraction is warranted. 

Comparison of the cavity phase of hydraulic ram to underwater ex- 
plosions provides valid relationships between cavityradius, energy and 


period of oscillation. 


B. PROPOSED RESEARCH 

HIGH SPEED FRAMING CAMERA 

A high speed framing camera capable of 200,000 frames per second 
has been obtained. When incorporated intothe enidtine test set-up, a 
complete time history of the shock front and projectile penetration will 
be available. 

STRAIN GAGES 

The existing entrance wall has been redesigned and fitted with several 
strain gages. Proposed research with this entrance wall should yield in- 
formation On wall deflections and wall bendingmoments caused by pro- 


jectile penetration. 
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APPENDIX A 


SAMPLE CALCULATIONS FOR THE GRID DISTORTION ANALYSIS 


fa fae Fa SES = Sas este line used for data 


nn te. ee line through 
center 





Fig. 12 Enlarged Grid Distortion 
Steps for obtaining a solution: 
1. evaluate equation (17) to find 2) 
2. evaluate equations (8) and (9) to find oe and Ap 
3. evaluate equation (12) to get & 
/ if 
4, evaluate equations (10) and (13) for © and 8 
5. evaluate equation (14) to get 4 
6. evaluate equation (15) to get Ur 
7, a equation (16) for Gh; 
e qu L 
8. evaluate equation (11) for - 
9. draw a curve of B vs r* using equation (18) 


10. using number pairs (B,r*) from step 9, evaluate equation (19) 


for the right hand side. Iterate until correct (B,r*) are found. 
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From Figs. ll and 12 


Qu 
il 


einen 
L+ = 19.8 inches 
L.,=9 inches 
“Yo = .35 inches 
fr = 3" 


D = 400 kg/cm? 
* A = (0.1) (180) = 0.29° = .005 radians 
(19. 8) 


From equation (17) 
0.35 = 35m(6,- 2) - {ye syaes (\-%)| TON 
Assume (®\-%) is small 
vy 0.352 3(0,-a) - (9-3)X 
3B, = 235+ IX 
or @, = 0.35 + 3% = .1318 radians 
S 


try 6, = 155° 


0- oO = 155° - 0.19= 1.1,- 
then 


0.BS= Este LL-[9- 3005 1,26] (tan 0.20) = 3486 


from iteration 
ee TSal ® 


. Sua b= Surb, (ae) = sm, (*\ 


t 
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Using p = 400 kg/om2 atr' = 3 inches 


e= IS” p = 9.00¢ 
W4= WR, ( [4 .004) = 1.39 (1,004) 


Ro [4353 
rR,. 133 ~ 06,9960 
Na 1,33255 

then 


sun b= bees aan) (.9960) 


Saby= -IBIL 

Pe b,= 1.2193° 

Ab= 64.-0,= 1.57 - 7.5133 
6b= 0,0517° 


from equation (12) 
= B- 4 © 

z= .19- ,051 = 0.0933" 

> 

from equation (10) 
@ = 90+ +4 6-8,= 904 0433 + .05(7-1,5] 
6'= 81.515" 

from equation (13) 

§" = 90+2) + 46- 6,= 94.515" 

from equation (14) 
Hi= T'SUACR-&)-Yy= Bsue (1,57- 118) — 035 
Y= 0,03" 
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from equation (15) 
Ve o[su1.41- gue (1,54- 19) 
Y= 0.015" 


from equation (16) 


$= 130 + .0>-1.57 


b/ = 111.4411 


from equation (11) 


= os AS Slot = -O>l] 





rele 
An approximation for B is 
\* = Yo +4) 
oe = 40+ = 0.25+49.03 _. 1167 
r » 


from step (9) 





B= 1.3353) 1308 .IU61) 
CILEINCI- IL6T) 


B=— 0.04947" 
from the assumed distribution for n 


ate ha + aay 
r' 


Se | oe are' + 0. OFO47 ( \-, 1167) when re- ais 
n= V Dalb> 
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then 


C= AN = WW 
Yr Vy 


€E= 1,3135- (,53 
[F55 
E= 0.0364, 


but 
6= IO AP 


". P= 107 (0.03646) 





Pp = 414.72 ey ene (calculated for a radius ratio . of Te 0 (LCT ) 
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